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Satellite system concepts that employ multiple spot beams with 
frequency reuse increase in complexity as the transmission capacity 
is increased and efficiently allocated throughout a service region 
much wider than the width of a single spot beam. A fundamental 
question is the ultimate capacity limit of a satellite system constrained 
in radiated power, bandwidth, and the aperture dimensions of the 
satellite antenna. Because the satellite system integrates multiple 
users, the capacity limit is expressed as a region over the space of 
transmission rates &\, •••, 3#g into the G ground terminals. Any 
particular set of rates contained in this region can be transmitted 
with arbitrarily high bit error-rate performance. In this paper, we 
find inner and outer bounds on the capacity region and show them to 
be dependent upon the geographical distances between ground sta- 
tions and the aperture dimensions of the satellite, since these control 
the effectiveness of frequency reuse. The outer bound is all-inclusive 
and no time-division multiple access, frequency-division multiple 
access, spread- spectrum, antenna-pattern shaping, or multiple scan- 
ning-beam concept can do better. Inner bounds are found by construc- 
tion. The bounds provide a yardstick against which the performance 
of any system can be compared. Detailed results of the effectiveness 
of frequency reuse are presented for the important case of two closely 
spaced ground stations. 

I. INTRODUCTION 

Communication satellites employing large antennas capable of pro- 
viding multiple, nonoverlapping, high-gain spot beams have attracted 

1089 



considerable attention in recent years because the allocated spectral 
bands can be reused in the various spot beams. 1 " 17 For a fixed level of 
available RF power, such satellites exhibit a much higher throughput 
capability relative to area coverage satellites because the higher an- 
tenna gain reduces the RF power required per spot beam to maintain 
the same effective radiated power (and therefore the same communi- 
cation performance) on the earth; the excess power can then be applied 
to the remaining spot beams, thereby increasing the capacity of the 
satellite. For widely separated spot beams, each such beam can be 
considered as providing an independent, full-bandwidth channel be- 
tween the satellite and the earth. A satellite offering N independent 
spot beams thereby provides a factor of N higher capacity than does 
an area coverage satellite. In fact, the total RF power required by the 
multiple spot beam satellite to provide this increased capacity is often 
less than that needed by the lower capacity area coverage system 
because the number of spot beams provided is usually smaller than 
the ratio of the spot beam to area coverage antenna gain. 

In this paper, we first briefly review several spot-beam, frequency- 
reuse concepts proposed to date. These include the familiar satellite- 
switched, time-division multiple access (ss-tdma) approach with fixed 
spot beams, 1 " 9 integration of the basic ss-tdma structure with a rapidly 
scannable spot beam to provide service over those regions not covered 
by a fixed spot beam, 10 " 15 and multiple scanning beam approaches 
which efficiently accommodate nonuniform service demand across 
some wide service region. 1417 This background material appears in 
Section II and is optional reading for those already familiar with these 
concepts. 

A fundamental question unanswered by these approaches concerns 
the ultimate capacity limit of frequency reuse communication satellite 
networks characterized by multiple ground stations arbitrarily located 
over some wide service area and by constraints upon the space platform 
radiated power, the available bandwidth, and the physical dimensions 
of the satellite's antenna aperture. This question is addressed from an 
information theoretic point of view in Section III, and inner and outer 
bounds on the capacity limit are derived. Because the satellite inte- 
grates multiple users with nonuniform traffic demands, the capacity 
limit is expressed as a region over the space of allowable transmission 
rates 0t\, • • • , 01 g such that independent information can be reliably 
transmitted at rate 0t x to ground station number 1, rate £? 2 to ground 
station number 2, ■ • • , and rate 01g to ground station number G. The 
G-dimensional region of acceptable rates &\, • • • , Mg defines the 
capacity region of the multiuser satellite system. For any traffic 
demand within this region, the power, bandwidth, and antenna aper- 
ture resources of the satellite can always be apportioned to satisfy the 
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demand. The derived bounds thereby provide a yardstick against 
which the capabilities of any system concept can be compared and also 
suggest areas in systems, coding, antennas, and devices where addi- 
tional research might greatly enhance the capabilities of future satellite 
systems. It is noted that the outer bound is all-inclusive; that is, no 
tdma, fdma, spread-spectrum, antenna-pattern shaping, or scanning 
spot-beam concept can provide a capacity region exceeding this bound. 
Attention is restricted to the down-link, although similar bounds can 
be derived for the up-link as well. We note that the Fourier transform 
relationship between the antenna aperture distribution and the far- 
field radiation pattern implies that the geographical distances between 
ground stations affects the capacity region. 

To find an outer bound on the capacity region, cooperation among 
the information sources is permitted at the satellite, and cooperation 
among ground stations is also permitted. Subject to radiated power, 
bandwidth, and satellite antenna aperture constraints, the total capac- 
ity is found for all combinations of G, G-l, • • • , 2, 1 ground stations in 
the network. For example, if there are two ground stations, the capac- 
ities into each individually are found, along with the combination of 
both treated as a single receiver. Call these Ci, C2, and C\,%. Then the 
capacity region outer bound is the intersection of the regions &\<C\, 

0t 2 < C 2 ,&\ +@2< Ci,2. 

By definition, any realizable approach for which the capacity region 
can be found provides an achievable inner bound to the capacity 
region. Transmission via independent beams in the absence of both 
cooperation among message sources at the satellite and cooperation 
among receivers provides one such bound. Each link from the satellite 
to a ground station is corrupted by Gaussian noise and possibly 
cochannel interference from neighboring beams. By varying the power 
division among beams, an inner bound is swept out. Time-sharing of 
several movable beams provides another inner bound. Tight inner 
bounds can often be obtained through judicious selection of both the 
number of beams chosen and their associated aperture distributions. 
It is noted that construction of the inner bound also provides a method 
for achieving the inner bound. 

Thus, in Section III we extend the concept of system capacity to 
multiuser satellite systems by extending the usual power and band- 
width constraints for a single-user channel to include an aperture-size 
constraint as well. The approach differs from earlier multiple-user 
analysis 18-21 in that here we are free to select the number of beams, 
their radiation patterns, and the apportionment of power among beams 
to maximize capacity. We note that these bounding techniques can 
readily be applied to a more restrictive satellite antenna wherein 
arbitrary aperture distributions subject only to aperture extent limi- 
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tations are replaced by other restraints. (For example, one might 
require that the aperture contain an iV-element, phased-array antenna 
of fixed element spacing and fixed elemental radiation pattern.) 

In Section IV, we invoke the theory of discrete prolate spheroidal 
sequences to numerically evaluate capacity region outer bounds for 
systems containing two closely spaced ground stations. This important 
case is examined in detail because, for two closely spaced ground 
stations, it has generally been assumed that time-sharing a single 
beam, rather than frequency reuse, must be employed to serve both. 
We can now, for the first time, investigate this hypothesis. Although 
analytical results can readily be obtained for systems containing three 
or more ground stations, such results cannot easily be displayed 
graphically (three or more dimensions are needed to portray the 
capacity region), and provide little physical insight into the utility of 
frequency reuse. Consequently, attention is restricted to networks 
containing two closely spaced ground stations. A typical result is as 
follows. Suppose the two ground stations are located at the — 1.3-dB 
contours (one on each side of center) of the narrowest beam which can 
be produced by the satellite. Suppose further that the maximum 
available carrier-to-noise ratio is 15 dB; that is, we could obtain a 15- 
dB carrier-to-noise ratio at one ground station if all available power is 
given to one beam with maximum gain in the direction of that ground 
station. Finally, both ground stations require the same transmission 
rate ^2 = 3.5 b/s/Hz. Then two simple pattern-shaped beams are found, 
one serving each of the two ground stations and each reusing the total 
available bandwidth, such that the power required is within 1.5 dB of 
the outer bound for that particular ground-station separation. Thus, 
the most complicated coding and radiation scheme possible could save, 
at most, 1.5 dB in required total power to achieve the rate ^2 provided 
by these two beams. Moreover, time-sharing a single scannable beam 
between the two ground stations would require 6 dB more total power 
to achieve the same capacity. 

In Section V, we qualitatively explore additional implications of the 
satellite system capacity region. Addressed here are the effects of rain 
fading, resource sharing to combat rain fading, and point-to-multipoint 
or broadcast satellite considerations. 

II. BACKGROUND 

The earliest attempts at employing multiple spot-beam, frequency- 
reuse concepts involved the formation of a number of up-link and 
down-link spot beams permanently directed to a corresponding num- 
ber of geographically isolated service regions which might appear as 
shown in Fig. la. 1 " 9 To interconnect the various spot-beam regions, an 
on-board time-division satellite switch is employed as shown in Fig. 
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Fig. 1 — Fixed-beam satellite switched time-division multiple access, (a) Typical fixed- 
beam service regions, (b) Time-division satellite interconnecting switch, (c) Typical 
interconnection sequences. 



lb; a time sequence of the various interconnections made by the 
satellite switch might then appear as shown in Fig. lc. Within each 
interconnection time slot, the ground stations of the network com- 
municate via time-division multiple access (tdma). Digital modulation 
consistent with tdma is employed. In this manner, the capacity of the 
satellite is time-shared among users. 
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To maintain cochannel interference among beams at an acceptably 
low level, the fixed spot-beam service regions must be spatially sepa- 
rated by perhaps two or more 3-dB beamwidths of the spot beam 
radiation pattern; closer spacing requires that a greater level of RF 
power be made available for each spot beam to overcome the degra- 
dation resulting from interbeam interference. Hence, the service re- 
gions are often separated by a much larger area within which no user 
can access the satellite (blackout region). 

Several approaches intended to reduce or eliminate the blackout 
region have been proposed. One technique consists of splitting the 
spectral band into several subbands and assigning these to spot-beam 
regions such that reuses of the same subband are separated by a 
distance sufficient to satisfy the interference criteria. 16 Additional spot 
beams are then provided using different subbands to cover the blackout 
region. Unfortunately, the demand for service can exhibit wide varia- 
tions among the spot-beam regions, and the band-splitting approach 
limits the available capacity into any one spot-beam service region to 
that provided by one subband. 

Another approach employs a pair of steerable spot beams (one for 
the up- and one for the down-link) in conjunction with fixed spot 
beams to provide service to the blackout region. 10 At each point in 
time, the up-link of the scanning beam can receive messages from a 
small spot-beam region contained within the blackout region, and the 
down-link can deliver messages to another small region within the 
blackout region. The scanning beam antennas appear as additional 
ports on the satellite switch. 1517 Then, by scanning the up- and down- 
links in such a manner that the amount of time dwelled at any one 
location is proportional to the traffic for that region, the traffic of the 
entire blackout region is accommodated by the single scanning spot- 
beam antenna. Thus, the scanning-beam concept combines the wide- 
access capability of an area coverage system with the power saving 
afforded by a spot-beam system. 

A drawback of the scanning plus multiple spot-beam concept arises 
from the fact that the demand for service from all regions (with the 
former blackout region now being treated as one large region) are 
dissimilar. Each region can now be served by the total available 
spectrum, but not all regions may require this capacity, implying waste 
of satellite resources. The total RF power available at the satellite can, 
however, be apportioned among the service regions on the basis of 
capacity requirements, thereby improving utilization efficiency some- 
what. 

Yet another approach involves the use of multiple scanning 
beams. 1417 Here, each antenna port is capable of providing a narrow 
spot beam that can be scanned anywhere over the total service region. 
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With this approach, it is possible to tailor the dwell time of a spot 
beam at any spot-beam region in proportion to the demand for service 
from that region. This concept is quite efficient in that a limited 
number of transponders are shared among a far larger number of spot- 
beam regions, thereby affording the power saving associated with 
highly directive spot-beam antennas while providing service over a 
wide region characterized by nonuniform traffic. The required number 
of spot beams is readily established by requiring that the total offered 
traffic (measured in time-slot units) equal the total number of available 
time slots. However, by employing a larger number of beams and 
transponders, each using the full bandwidth, one can envision the use 
of channel coding techniques within each beam to save power by using 
the available excess bandwidth. The number of scanning beams which 
can be provided before interbeam interference becomes a problem is 
still an open question, however. Because the dwell time at any one 
location is dependent upon the traffic from that location, this number 
must be dependent upon the terrestrial traffic distribution. 

Thus we see that various approaches have been proposed to intro- 
duce spot-beam frequency-reuse techniques for satellite communica- 
tions. These represent varying degrees of refinement of the basic 
approach in an attempt to more efficiently allocate the power and 
bandwidth resources of the satellite such that a uniform grade of 
service is maintained over a service area much larger than the beam- 
width of a single spot beam. The traffic patterns might be highly 
nonuniform over the total service area. 



III. DERIVATION OF THE CAPACITY REGION BOUNDS 

The satellite system to be studied is shown in Fig. 2. Here, the 
satellite is presented with G independent messages Mi through Mg to 
be communicated, respectively, to geographically separated identical 
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Fig. 2 — Satellite system model. 
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ground stations 1 through G. The receiving gain of each ground station 
is normalized to unity. These messages are typically produced at the 
satellite via onboard detection of the up-links. Thus, message source 
Mi, \<j-< G, represents the composite of all up-link communication 
from the various ground stations intended for reception at ground 
station number/ It is assumed that the various up-link messages are 
detected at the satellite error-free. This is consistent with up-link 
transmission at rates lower than those permitted by the up-link capac- 
ity region. Since power constraints are usually not as tight on the up- 
link as on the down-link, we expect the satellite system capacity to be 
dominated by the down-link capacity region, and bounds are found 
only for the down-link; similar bounds can be established for the up- 
link, however. Message source Mj produces information at a rate of 
Stj b/s, j = 1, • • • , G. The satellite is constrained to a total radiated 
power of P watts and an RF radiation bandwidth of W Hz. The 
antenna aperture distribution is identically zero outside some two- 
dimensional region in the plane of the aperture. 

3.1 Outer bound 

We seek first an outer bound on the collection of rate vectors^? = 
\0t\t • • • , ^?g] for which reliable communication from the satellite to 
the ground stations is possible. Within this region, the power, band- 
width, and aperture dimension resources of the satellite can be allo- 
cated to achieve error-free communication at rates ^. 

Let message source no. i produce the binary information stream a,, 
i = 1, • • • , G. The collection of information streams can then be 
represented by the matrix g = [fli, • • • , flo]. For a particular message 
matrix & the aperture distribution at the satellite can be represented 
as 

sf(x', t, a) = Re[s?(x', t, a)e J "°'l (1) 

where wo is the carrier frequency and j/ is the complex baseband 
aperture distribution. In the above, we have, for simplicity, assumed a 
one-dimensional antenna; the approach can be generalized to a planar 
antenna, however. The far-field radiation pattern S(x, t, a) is obtained 
via a Fourier transform of the aperture distribution. For baseband 
bandwidth small compared to the carrier frequency w , 



S(x, t,a)=Re 



— j j£(x' t t, & )e ikxx ' /Zo e Juiot dx' 

2.TT J 



(2) 



where k = w /c, c is the speed of light, and z is the distance between 
the plane of the satellite and the plane of the earth. Equation (2) 
applies in the far field, i.e., z » aperture extent. 
Let 
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A(x t, a) <& j2/(x', t, a), (3) 

where (<=>) implies the Fourier transform. Then the farfield radiation 
pattern is: 



S(x, t, a) = A r (x, t, a)cos co t + A t (x, t, a)sin u t, 



(4) 



where A r and Ai are, respectively, the real and imaginary components 
of A 

We see that the far-field radiation contains in-phase and quadrature 
components. Again, because the baseband bandwidth is small com- 
pared to the carrier frequency, these two components are orthogonal. 
Thus, the in-phase and quadrature components can provide two in- 
dependent, identical capacity channels. Attention is therefore re- 
stricted to one channel, say, the in-phase channel; results for the 
quadrature channel are identical. 

The ground stations are located at Xi, • • • , x G (again, for simplicity, 
we assume that ground stations are colinear; the approach can be 
extended to planar constellations, however). The ground stations are 
allowed to communicate via a noiseless terrestrial link to a central 
processor as shown in Fig. 3. At the ith ground station, we receive 

<W(Xi, t, a) = A r (Xi, t, a)cos u t + m(t)cos u t, (5) 

where m(t) is white, normal, and of spectral density N /2. We see, 
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Fig. 3 — Processing model for capacity region outer bound. 
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then, that the central processor is presented with the G-dimensional 
received vector 



M = 



A r (xi, t, a) + ni(t) 
A r (x 2 , t,a) +n 2 {t) 

_A r {x G , t, a) + n G {t)_ 



cos Uot, 



(6) 



where m, »y are independent, i f*j. 

Let us express A r (xi, t, a) as an expansion in some complete ortho- 
normal basis {u(t)}\ 

AAxi, t,a)=X bj( Xi , a)uj(t). (7) 

j 

Then, the part of the received vector relevant to detection is 



bj(xi, a) 

1 Ibjixo, a) J 



uj(t) + % 

J Lna/J 



Uj(t) \ COS CJo^, 



where 



mj = rii(t)uj(t) dt, i = 1, 



G. 



Clearly, 



E{nijn k , P } = — 8i, k 8j, p . 



(8) 



(9) 



(10) 



Let vectors Ji, • • • , Jg be an orthonormal basis spanning G-dimen- 
sional Euclidean vector space. Then, for a particular set of ground 
station locations Xi, • ■ • , x G , and for every message matrix a, we can 
expand Z in terms of the orthonormal basis {J) : 

Z = { S I Kij(a)uj{M + X S -r\ijUj{t\£i [cos w *, (11) 



i ,=\ 



where 



and 



Vti = 



•li 



Ln-GjJ 



(12a) 



(12b) 



In (11), the expansion coefficients K,j{a) are dependent upon the 
particular basis {€} employed. 
Thus, we see that since reception of the radiated signal can occur at 
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no more than G points, the total channel from the satellite to the 
ground stations consists of G independent parallel channels, each of 
bandwidth W. Viewed from the perspective of the G ground terminal 
receivers, the satellite is, in effect, radiating G-dimensional vectors 
rather than continuous functions of space since reception can only 
occur at the G points where receivers are located. We conclude that 
the most general baseband radiation to maximize the capacity of the 
G parallel satellite channels is of the form 



a(t, x,a)=n LAa)Ux)uj(t), 



(13) 



; '=i 



and the signal portion of the received vector is of the form 






In (14), the vectors »£, = 



Uxg) 



LMxg) J 

and ^j = 



Uj(t). 

Uxa) 



(14) 
are orthogonal, 



i ¥* j. For convenience, we normalize the energy of the underlying 
radiation patterns to unity, i.e., / | ifc(x) \ 2 dx = 1, i = 1, • • - , G. The 
functions \l>i(x) and ty(x) are not necessarily orthogonal, however, over 
the real line, i j^ j. 

We note from (14) that some radiation patterns will deliver more 
radiated power to the G ground stations than others. The problem, 
then, is to find G radiation patterns ^(x), producing G orthogonal 
receiving vectors i£,-, which maximize the total capacity, subject to the 
constraint imposed by the finite radiation aperture of the satellite 
antenna. This maximization is performed under the additional con- 
straint that the total radiated power is bounded. Since the G receiving 
vectors are orthogonal, selecting their inputs to be independent maxi- 
mizes the capacity, i.e., 



E{$utk*}-E{£i3}8i, k . 



(15) 



Let the total power radiated by the satellite be constrained to P 
watts. Thus, for each a, 



i 2 (t, x, a) dt dx < P, 



(16) 



where T is the temporal extent of the transmission. Since we are 
interested in obtaining an outer bound to the capacity region, any 
bound derived from a constraint which contains (16) as a special case 
must be an outer bound to the true capacity region. Rather than 
insisting that (16) hold for each a, let us impose the constraint that, 
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averaged over all information matrices, the radiated power must not 
exceed P, i.e., 

^f I (At,x,a)dtdx) <P. (17) 

Substituting (13) and invoking the orthonormality of the set {u(t)}, 
the normality of the set {\p(x)}, and the independence of the encoded 
symbols {£}, we obtain: 

£ 2 <fi/> = P. < 18 > 

i «=i 

For each i, we can interpret £/ < £?, y > to be the power in the ith 
mode. For convenience, let 

S <fl/> = «iP < 19 > 

y 

| a, = 1. (20) 

Each mode can radiate both in-phrase and quadrature components 
of the received process. Thus, the ith mode provides an additive white 
Gaussian noise (awgn) channel of bandwidth W and power on P. 
Invoking the classical formula for the awgn channel, the total capacity 
for the G parallel channels is seen to be 

r mP 1 4>h] 

c *-l w M 1+ ^Br-\> (21) 



where i//,,* is the Ath component of the vector 

_^,(*g)J 



(22) 



and \pi(x) is the normalized radiation pattern corresponding to the ith 
mode. In writing (21), we have assumed ground-station receiving 
apertures normalized to unity. This normalization is permissible since 
we will soon rewrite (21) in terms of the available carrier-to-noise ratio 
(cnr) at a receiving terminal. We note that, for all i, k, the vectors ±i 
and ±k are orthogonal, i.e., 

± i .± k =\± l \%k, (23) 

but that the functions ty(x) and \p k (x) need not be orthogonal. 

For each G-member set of functions {^(x)} and for each power 
allocation oti, "",ac, the total capacity is given by (21). We note 
that this can be maximized by allocating power according to the 
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classical parallel channel water-fill problem. 22 Here, we consider that 
G Gaussian parallel channels are available, of noise powers 

, AW 



o, = 






1=1, 



G. 



(24) 



ti« 



The total available power is P, and we must select a\ , • • • , olg satisfying 
the constraint (20) such that the capacity is maximized. The solution 



is 



a 2 + aiP = B for of < B 
a, = for a 2 > B, 



(25) 
(26) 



and B is chosen such that (20) is satisfied. This solution has the 
geometrical interpretation shown in Fig. 4. Available is a quantity of 
water P to be placed in a vessel with a bottom consisting of J levels as 
shown (J = 6 in Fig. 4). The height of the ith level is the noise power 
of. Then, the amount of power given to the ith channel is equal to the 
height of water above the ith level. For any channels exhibiting too 
high a noise power for the available signal power, the greatest capacity 
is achieved by allocating no power to those channels. 

The G normalized functions xfriix), • • • , \^g(x) were chosen such that 
(23) is satisfied for all i, k and are also aperture-limited, that is, the 
Fourier transforms of these functions vanish outside the finite dimen- 
sions of the antenna aperture but are otherwise arbitrary. To maximize 




S Pi = p 

i»l 



Fig. 4 — Parallel channel water-fill analogy. 
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the capacity of the G parallel channels, we must select those modes or 
functions \pi(x), • ■ • , \p G (x) satisfying the constraints and must further 
allocate power among the modes such that (21) is maximized. It is 
convenient to first place (21) in a form involving available cnr at a 
receiving terminal. 

Consider an antenna aperture limited in extent to [-A, A]. The 
maximum gain available from this antenna in the direction of a ground 
station located at x = arises for the normalized aperture function: 
, i!)( . , |= , , r ,-A<x'<A (2?) 



i^jA, 
c ' |0, else 



For this function, 



*(*)« 



elsewhere. 



[A sin kAx/zo 
•n kAx/zo 



(28) 



Thus, if all power is given to this mode, the power available at a 
ground station located at x = is Pg max , where 



_A 

IT 

For this case, the carrier-to-noise ratio is simply 

x gmax 



P = 



NqW 



(29) 



(30) 



Using the relationship (30), the total capacity for the G ground 
terminal network is given by 



Ct = supremum < £ W 1°S 2 



where 



1 + p<Xi 



G 

I 



I 4b 



gr 



tf(kx/z ) d(kx/z ) = 1, i = 1, • • • , G 



(32) 



(33) 



^(kx/z )e- Jxx ' /z °d(kx/zo) = 0, | x' | > A (34) 



L<M*g)J 

G 



(35) 

(36) 
(37) 



_A 

gmnx • 



(38) 
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Clearly, the location of the ground stations Xi , • • • , xg and the antenna 
aperture 2A affect the maximum capacity. We note from (36) that a, 
and 07 are included in the orthogonality constraint because the radia- 
tion vector of any mode receiving zero power (i.e., a = for that mode) 
need not satisfy any orthogonality conditions. 

Condition (34) implies that, to maximize (32), each \pi(x) can be 
expressed as some linear combination of G basis functions, each itself 
a linear combination of the G sampling functions [sin kA(x — xA/zo]/ 
kA(x — x/)/zo] € — 1, • • • , G, with the property that the £th basis 
function equals unity at x = Xh and equals zero at x<, € V k. 23 The 
coefficients in the expansion of each \pi(x) are then chosen, subject to 
(36), such that (32) is maximized. When the ground stations are 
uniformly spaced, the theory of discrete prolate spheroidal se- 
quences 24,25 can be applied to more directly yield the capacity. 

Having found the total capacity possible for a G-terminal network 
with power, bandwidth, and aperture dimension constraints and with 
noiseless channels connecting each ground station to a central proc- 
essor unit, we conclude that the sum of the information rates into the 
G-terminals cannot exceed Ct, i.e., 

^1+^2+ • • • +0t G <C^. (39) 

The capacity region must therefore be totally contained within the 
region (39). 

Proceeding in an identical manner, we can find the maximum 
capacity possible into a network containing any G-\ of the original 
ground stations, assuming that all of the available satellite resources 
are available for communicating to the G-l terminals chosen. In this 
manner, we find G capacities Cr7m, rn = 1, •••, G, one capacity for 
each choice of G-l out of G participating ground stations. We conclude 
that the capacity region must be totally contained within the region. 

0t x + • • • + 9tG-\ < C% 1 (40) 

0t\ + ••• +^ G -2+^G<C^ 1 (41) 



+ m • . • + 0t G < Cf7d-i (42) 



#2 + #3 + • • • + #g=£ Cf; G \ (43) 

Continuing, we find the maximum capacity for all combinations of 
G-2, G-3, • • • , 2, 1 ground stations in the network, and write an 
inequality for the sum of the rates for each combination. An outer 
bound on the capacity region for the G-terminal network is then the 
intersection of the regions corresponding to each inequality. For ex- 
ample, if G = 2, we can find three regions, 
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h + m 2 < Ct 

St\ < C T ,i 

^?2< Cfo. 



(44) 
(45) 
(46) 



The intersection of these three regions, shown in Fig. 5, is an outer 
bound on the capacity region of the two-terminal network. 

We note that the outer bound (32) can be applied to other aperture 
constraints more restrictive than the size constraint (34). For example, 
one might require that each member of the set {^(x)} be realized by 
a phased-array antenna with a specified-element radiation pattern and 
a fixed-element separation. Then each member of the set {\p(x)} would 
be represented as a linear combination of the radiation patterns of the 
elements, and Ct would be represented as the supremum over the 
coefficients in each such linear expansion and the set {a}. The coeffi- 
cients must then satisfy a normalization relationship such that (33) is 
satisfied. Alternatively, the antenna might be specified as being of 
offset Cassegrainian design with fixed-feed arrangement, and each 
member of the set ty(x)} would be represented as a linear combination 
of the radiation from each feed. Other constraints are, of course, also 
possible. For any other set of constraints, however, the resulting 
capacity region would be smaller than that obtained via (32) through 
(38), i.e., the capacity region obtained via (32) through (38) contains, 
as a subregion, the capacity region obtained by a constraint more 
restrictive than a size constraint. 




Fig. 5 — Typical outer bound on the capacity region for two ground stations. 
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3.2 Inner bound 

We now turn our attention to obtaining a lower bound on the 
capacity region of the G-terminal network. To find a lower bound, we 
remove the central processing unit of Fig. 3 as well as the noiseless 
terrestrial interconnecting links. Any method for then communicating 
the informational contents of message sources M\ through Mg to the 
intended destinations provides an inner bound to the capacity region 
if the capacity region for that method can be calculated. The tightness 
of the bound varies with the particular method chosen. 

One method for communicating between the satellite and the ground 
stations consists of time-sharing a single channel and antenna beam 
among all destinations. The antenna can form a single beam, directed 
anywhere, of maximum gain g mi , x ; when focused at one particular 
receiving terminal, the available carrier-to-noise ratio is p. The single 
channel capacity for this transmission scheme is then 

C, = Wlog[l + p]. (47) 

By varying the amount of time that this single beam dwells at each 
ground station in proportion to the traffic demand for that station, we 
obtain the following lower bound for the capacity region: 

:^, + . • • + M G < C, = W\og[l + p]. (48) 

Another inner bound can be found by having the satellite form G 
beams, each of maximum gain in the direction of one ground station. 
The ith beam is given a,P watts of power, Yli"=\ a « = 1. To transmit at 
capacity, the encoded input symbols to each beam are independently 
drawn from a Gaussian distribution of variance a,P, i = 1, • • •, G. The 
spectrum is totally reused in the various beams, and the channel inputs 
to all beams are independent. 

The pth ground station observes both the desired signal as well as 
independent additive Gaussian interference from all other beams; the 
strengths of these interferers are dependent upon the locations of the 
ground stations. Let k p ., be the power from the 1'th interferer observed 
at the pth ground station relative to the desired signal at the pth 
ground station. For maximum gain beams, 

sin{kA(x p - Xi)/z ) 

{kA(x p -Xi)/z„} 

i= 1, .--, G; p = 1, .-., G. (49) 

The rate at which information can be reliably communicated into the 
pth ground station is then bounded by 

a p p 



m p < W\o% 



1 + 



a 

1 + P H (XiKpj 






p=l, ...,G. (50) 
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For a fixed allocation of power to the various beams, the point ^ 
obtained via (50) with equality provides an achievable rate vector and 
is therefore a point on an inner bound to the capacity region. By 
varying the power allocation among beams, an inner bound on the 
capacity region is swept out. 

Yet another inner bound is obtained by again forming G independent 
beams, each of maximum gain in the direction of one ground station 
subject to the constraint that no interference be radiated in the 
direction of any other ground station. The encoded channel inputs to 
each beam are selected from independent Gaussian distributions. 
Thus, G independent, noninterfering channels are created. Again, the 
noninterfering constraint implies that the maximum gain for each 
beam formed is dependent upon ground- station locations. Non-inter- 
fering beams are found as follows. 

Consider the pth ground station. Let Y p {x') be the aperture distri- 
bution for the pth beam. Then the radiated field for this distribution 



is 



u p (x) = — Y p (x')e> xx dx', 



(51) 



where, for convenience, we have normalized the aperture dimension A 
and the scale factor k/z a both to unity. Now we require that 

u p (xj)=0, j¥>p. (52) 

Using Lagrange multipliers aj,j ¥* p, we wish to maximize 



Up(x p )\ 2 = 



£ |_ *'<*'> 



G 

j-i 
Mp 



dx' 



(53) 



Invoking the Schwartz inequality, we find the maximizing aperture 
distribution to be 

i> p (x') = B e- jx " x ' - I aje- jx J x ' , (54) 

Mp 

where the constants [a,, j ¥* p) must be selected to satisfy the 
constraints (52) and B is selected such that f \ u p (x) \ 2 dx = 1. From the 
constraint equations (52), 

•l r g 



B 

u p {xj) =0= — 

Z7T 



e -J*px- _ £ a k e- Jx k x ' 
k=i 

k¥*p 



e JX ' x dx (55) 



stt\{xj - x p ) £ sinfo - x k ) , „ 
= 2, a k > y = l, •••,(j. 

Xj — X p k-l Xj Xk J ^k p 

k*p 



(56) 
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The G-l equations (56) can be solved for the unknown parameters a,, 
j ^ p, and the normalization factor B follows. The remaining G-l 
noninterfering beams are found similarly. For a given power allocation 
{a}, we then obtain 



?;< W\og 



1 + CLjP 



gj(Xj) 

gmax 



7 = 1, ---.G, 



(57) 



where gj(xj) < g max is the gain of they'th mode at the location of the 
y'th ground station. Again, by varying the power allocation, we sweep 
out a lower bound on the capacity region. 

The three inner bounds for G = 2 might appear as shown in Fig. 6. 
We note that, depending on the ground station locations, the multi- 
beam bounds may or may not exceed the time-share bound. Other 
inner bounds can also be found. Of course, none of these bounds can 
exceed the previously found outer bound. 

We note that, by construction, derivation of an inner bound yields 
a method whereby the inner bound can be attained since, when 
deriving these bounds, we did not require cooperating ground termi- 
nals. 

IV. DETAILED STUDY OF TWO GROUND TERMINAL SYSTEMS 

The bounding techniques introduced in the previous section are now 
applied to study the capacity region for two ground terminal systems 
subject to power, bandwidth, and antenna aperture dimensional con- 




V^\ BOUND FOR 

x . \ MULTIPLE MAXIMUM 

\ \ „_-GAIN BEAMS SUBJECT 
\ N ^" TO ZERO INTERFERENCE 

\ \ CONSTRAINT 

\ 
\ 
\ 

TIME-SHARE _-- Jr \ ^ 

bound" 



V 

\ 



BOUND FOR 
, -MULTIPLE MAXIMUM 
<*" GAIN BEAMS 

\ 

\ 

\ 
\ 



1Nl_ 



C, = Wlog 2 (1 +p) 



C, R, 



Fig. 6 — Three typical capacity region inner bounds for two ground stations, p is the 
resulting carrier-to-noise ratio at a receiving ground terminal if all RF power available 
at the satellite is contained in a single beam of maximum gain in the direction of that 
ground terminal. 
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straints. Solution to this problem amply illustrates the relationship 
between ground-station separation and frequency-reuse advantages, 
while providing results that are easily presented and evaluated. More- 
over, a study of the capacity region for two closely spaced ground 
stations provides, for the first time, a comparison of frequency-reuse 
techniques and simple time-sharing on the basis of available capacity. 
The analytical techniques used to evaluate the outer bound (32) are 
readily applied to study networks containing an arbitrary number of 
equally spaced ground terminals. A more general technique, described 
earlier in Section 3.1, can be applied to study systems in which the 
ground terminals are arbitrarily located. 

To study the outer bound for a two-ground terminal network, we 
invoke the theory of discrete prolate spheroidal sequences (dpss) 24,25 
to find the supremum of (32). This theory is reviewed in Appendix A. 
Using the theory and notation contained therein, we now consider the 
expression for capacity (32). The G capacity-maximizing normalized 
radiation patterns {i//,} are aperture-limited via (34) to [—A, A]. Invok- 
ing the sampling theorem, we can express \f>i(x) as 

\pi(x) =Xo }j winXo) ; ttt , (58) 

„=-«, tt\x — lULo) 

where 

X = ■%-, A > A. (59) 

A 

In (58), we have again normalized the scale factor k/z a to unity since, 
as will become apparent, the capacity region bounds are dependent 
upon the separation between ground stations relative to the null 
spacing of the narrowest beam which can be produced by the satellite 
aperture, rather than upon absolute distance. 
Since J$(x)dx= 1, 

i tf(nX )=^r. (60) 

Thus the normalization (33) implies the constraint (60) on the sample 
values of the radiation. 

Let the G ground stations be located at X , 2X , ■ • •, GX . The 
expression for capacity (32) depends only on the vectors [\pj(X ), 
\pj(2X ), ••• ^j(GXo)], j = 1, •••, G. Each of these vectors must 
therefore be linear combinations of the index-limited dpss, or else a 
new set of vectors satisfying (36) could be found with greater energy 
at the location of the ground stations. 

Restricting attention now to two ground stations, we can write 

fr (nX ) = avi (nX ) + b v 2 (nX ) (61) 
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t// 2 (nXo ) = cvx (nXo ) + dv% (nX a ) , (62) 

where vi# = [• • • pi,2(— X,), 1*1.2(0), Vi^{X ), 1*1,2(2X0), • • •] are propor- 
tional to the dpss for w = tt(A/A ), N = 2 (the proportionality 
constant is (1/X ) 1/2 ). From (60) and (95), we find 

a 2 + b 2 = 1 (63) 

c 2 + d 2 = 1. (64) 

Furthermore, from (96) 

£ tffciXo) = 2 [a 2 v\(nX ) + b 2 vl(nX )] = -[<*% + 6 2 A 2 ] (65) 

n— 1 n=l *o 

2 ^KnXo) = I [cM(nXo) + tf 2 i> 2 2 (/iX )] = i-[c 2 A, + d%]. (66) 

n-l n=l Ao 

Applying (36) and (96), 

a{l - a)[\iac + \ 2 bd] = 0. (67) 

Thus, applying (38), (32) reduces to finding a, b, c, d and a such that 

C9 = W\og[l + gi ap] + Wlog[l + g 2 (l - a)p] (68) 

is maximized, subject to (63), (64), and (67), where 

ft _«!X^6% m 

g2=£ !^% (70) 

In the above, D = A/A < 1 is the separation between ground station 
expressed as a fraction of the first null spacing of a single-mode 
radiation pattern possible with the aperture of extent 2A. Figure 7 is 
the radiation pattern of a single maximum-gain mode directed at one 
ground station. This radiation possesses the familiar (sin X)/X pattern. 
For D = 1, the two ground stations are located as shown. 

The eigenvalues Ai and A 2 of (92) are dependent only upon D: 

n sin ttD 
A,=D + (71) 

77 

^ sin ttD 
\ 2 = D . (72) 

7T 

For two ground stations, an outer bound for the capacity region is 
given by the intersection of the three regions 

®x + 0ti < Cf (73) 

^i,2 < Wlog[l + p]. (74) 
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RADIATION VOLTAGE 



MAXIMUM GAIN BEAM 

.,-— — sin x/x 




Fig. 7-*-Ground station locations defining the condition D = 1. By definition, D is the 
separation between ground stations expressed as a fraction of the beam center to first 

null of the narrowest beam that can be produced by the satellite antenna aperture. 



Figure 8 shows the outer bounds for select ground station separations 
D and select carrier-to-noise ratios p (the rates have been normalized 
to b/s/Hz). In Fig. 9, we present additional bounds for lower values of 
p. For sufficiently high p, these curves clearly show that, as the distance 
between the two ground stations is reduced, the capacity region 
becomes smaller. At smaller ground-station separation distances, fre- 
quency reuse becomes less effective since, to maintain orthogonality 
of the modal received vectors, a smaller fraction of the energy of the 
second mode is directed toward the two ground stations. The water- 
fill solution then implies that, as the ground stations approach one 
another, a greater fraction of the available power is allocated to one 
mode to maximize the capacity. 

For lower values of p, however, we see that the capacity region outer 
bound increases as D becomes smaller. This apparent anomaly is 
readily explained. For small p, there is insufficient power to excite 
more than one mode, and frequency reuse is of no advantage. As the 
distance between ground terminals is reduced, more power from this 
single mode is intercepted by the ground terminals, and the capacity 
region thereby increases. 

Figure 10 plots the mode radiation patterns for a carrier-to-noise 
ratio p = 9 dB and a separation D = 0.3. In this case, 90 percent of the 
available energy is allocated to one mode. As the carrier-to-noise ratio 
increases, the two radiation modes remain the same but the power 
allocation changes such that, at p = 21 dB, both modes receive the 
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same energy. Conversely, for values of p smaller than 9 dB, mode 2 
receives no power at all. 

Plotted in Fig. 11 as functions of separation D are the carrier-to- 
noise ratios below which the capacity is maximized by one-mode 
operation and above which the power allocation imbalance is less than 
10 percent; that is, both modes receive approximately equal power. As 
D approaches unity, the cnr for each approaches — oo. 

Figure 12 plots an inner bound to the capacity region for select 
values of p and D. This bound is a composite of (a) the time share 
bound and (b) the bound obtained by exciting two noninterfering 
modes. For certain (^1,^2) pairs, the greater inner bound is achieved 
by time-sharing; for other pairs, the opposite holds. For G = 2, we 
solve eq. (56) and normalize to find [see eq. (57)]: 

gl (X„ ) = g2 (2X„) = 1 - (SiD " D/ * D) \ (75) 
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Fig. 8 — Two ground-station, capacity-region outer bounds for select carrier-to-noise 
ratios p and select ground stations separations D. See captions for Figs. 6 and 7 for 
definitions of p and D. 
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Fig. 9 — Two ground-station, capacity-region outer bounds for select carrier-to-noise 
ratios p and select ground station separations D. See captions for Figs. 6 and 7 for 
definitions of p and D. 



In Fig. 13, we consider the special case where the two ground stations 
require equal capacity {0t\ = 0li = 0t). Two sets of curves are plotted. 
The first set shows, as a function of ground separation D, the difference 
in carrier-to-noise ratio between the upper bound and the noninterfer- 
ing dual mode lower bound to maintain select transmission rates 0t 
into each ground station. From this curve can be inferred the maximum 
power saving possible via the most complicated coding and radiation 
schemes to achieve a rate 0t, compared to interference cancellation. 
The second set of curves shows, as a function of separation D and for 
select values of ^2, the power difference between time-sharing and 
interference cancellation to maintain that value of 0t. For example, 
consider the case 0t = 3.5 b/s/Hz, D = 0.5. To reliably maintain this 
transmission rate into each ground station via interference cancellation 
requires a carrier-to-noise ratio of 15 dB. The upper bound shows that, 
to maintain this rate, the minimum required cnr is 13.5 dB. Thus, the 

1112 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1 980 



power saving possible by the most complicated coding and radiation 
scheme cannot exceed 1.5 dB. By contrast, for simple time-sharing, a 
CNR of 21 dB is needed. Thus, the method of interference cancellation 
achieves performance almost as good as the best theoretically possible 
and provides a power saving of 6 dB relative to time sharing. As shown 
in Fig. 14, for D = 0.5, both ground stations are well within the -3-dB 
contour of a single maximum-gain beam and therefore one might 
attempt to serve both either by the same fixed beam or by scanning a 
single beam between them; the results here show that it is possible to 
serve the two ground stations by two independent beams totally 
reusing the available spectrum and thereby save 6 dB in satellite 
power. 

V. EXTENSIONS 

The bounds on the capacity region for frequency reuse point-to- 
point communication satellites derived in Section II can be applied to 
other situations. We now discuss, qualitatively, two such extensions, 
confining attention exclusively to the outer bound. 




RELATIVE DISTANCE- 



Fig. 10— Two mode radiation patterns for p = 9 dB and D = 0.3, weighted by optimum 
power allocation. See captions for Figs. 6 and 7 for definitions of p and D. 
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Fig. 11 — Maximum CNR for single-mode operation and minimum CNR for equally 
excited dual mode operation vs ground-station separation. 



At frequencies above 10 GHz, satellite communications are strongly 
influenced by precipitation. 26 The major effect of precipitation at these 
bands is the attenuation or fading of signal strength by heavy rainfall. 
The bounds of Section II were derived under the assumption that 
signal is available at all ground stations unattenuated. Thus these 
bounds always apply at frequencies below 10 GHz and apply for clear 
air conditions at higher frequencies. The bounds, however, can readily 
be extended to account for rain-induced fading. 

Suppose that, due to rainfall, the signal voltage at they'th ground 
station is reduced to fyV, where V is the clear air value and < f3 < 1. 
Typically, at any given time, /?> = 1 at most ground stations. For a 
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particular set {/?,} , an upper bound to the capacity region for the G 
terminal network is then 

G 

g P a J 2 PUi*" 1 

I W\og 2 



G _ supremum 

T "{«},WW) 



J-\ 



1 + 



A-l 



gn 



where 



G / G 

}2.i. .1. _ ,21 v o2 |2 



otjOii I tffeifo = aj 2 W«l «/./ 



A=l 



A=l 



(76) 



(77) 



and conditions (33) through (35), (37), and (38) apply. Equation (77) is 
the orthogonality condition among modal received vectors, replacing 
(36), applying in the presence of rain fades. 

The result (76) has an interesting interpretation. Suppose that, for 
clear air conditions, the vector^? = {^?i, • • • , Mg\ at which information 




(BITS/SECOND/HERTZ) 

Fig. 12 — Two ground-station, capacity-region inner bounds for various carrier-to- 
noise ratios p and ground station separations D. See captions for Figs. 6 and 7 for 
definitions of p and D. 
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Fig. 13 — Difference in decibels between the achievable capacity for two maximum 
gain, noninterfering modes and (a) the capacity region outer bound and (b) the time- 
share inner bound vs ground-station separation distance for selected equal rates into 
both ground stations. See Fig. 7 for definition of D. 

is communicated to the ground stations is not on the boundary of the 
clear air capacity region. For various sets of fading conditions, corre- 
sponding to different attenuation factor sets {/?} , we calculated the 
degraded capacity region via (76). Then, provided i# is within the 
degraded capacity region for all fading conditions of interest, we can 
reallocate the satellite power, bandwidth, and antenna resources 
such that reliable communication is maintained at the rates PA_. For 
example, to maintain the desired service availability, it may be neces- 
sary that up to F simultaneous fades, each of fade depth k dB, be 
accommodated. There are G\/[F\ (G - F)\] ways that the F fades can 
occur at the G ground stations. For each way, we find the resulting 
capacity region. Then, provided ^ is within each capacity region, an 
allocation can be made to maintain service. For the typical case where 
F «: G and all ground stations require approximately the same infor- 
mation rate, the vector.^ can be close to the boundary of the clear air 
capacity region. These observations are the essence of resource sharing 
approaches to combat rain fading, discussed elsewhere, 27 in that by 
accepting a traffic flow slightly smaller than the available clear air 
capacity, we can reallocate satellite resources to maintain service 
during fade events. 
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The capacity region outer bound developed in Section II can be 
applied to point-to-multipoint or broadcast satellites as well as point- 
to-point systems. For such an extension, we consider that, at the 
satellite, the G independent message sources, each one of which is 
intended for reception at a single ground station, are replaced by 
2 G — 1 independent message sources divided into G groups containing 
G\/[k\(G — k)\] members, k = 1, • • •, G such that each source within 
the £th group generates information intended for reception at one 
combination of k ground stations out of the total of G ground stations. 
For example, there are three groups for G = 3 ground stations. Three 
sources M\, M-2, and Ms within the k = 1 group produce information 
at rates ^?i, .^?2, and 0t% intended for reception at ground stations 1, 2, 
and 3, respectively. Three sources M ii2 , M\, 3 , and M 2 , 3 within the k = 
2 group produce information at rates St\^ % ^1,3, and ^2.3 intended for 
reception at ground stations 1 and 2, 1 and 3, and 2 and 3, respectively. 
One source M\ ,2,3 within the k = 3 group produces information at the 
rate 0t 1,2,3 intended for reception at all three ground stations. When 
deriving the outer bound, it was assumed that ground stations can 
communicate via a noiseless terrestrial link to a central processor. 
Thus, using (32), we can bound the rates of the seven broadcast 
satellite sources for G = 3 ground stations as follows: 




^GROUND STATION 
LOCATIONS 



DISTANCE D- 



Fig. 14 — Ground-station separation for D = 0.5 relative to a single maximum gain- 
beam radiation pattern. See caption for Fig. 7 for a formal definition of D. 
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Sty + #2 + #3 + #1.2 + #1.3 + #2.3 + #1.2.3 < Cf (78) 

#, + #2 + #1.2 + #1.3 + #2.3 + #1.2.3 < Cf(l % 2) (79) 

#, + #3 + #1.2 + #1.3 + #2.3 + #1.2.3 < Cf (1, 3) (80) 

#2 + #3 + #1.2 + #1.3 + #2,3 + #1.2.3 < C? (2, 3) (81) 

#! + #1.2 + #1,3 + #1.2,3 < CV (1) (82) 

#2 + #1.2 + #2.3 + #1,2.3 < CP (2) (83) 

#3 + #1.3 + #2.3 + #1.2.3 < CP (3). (84) 

where CP(i,j) is the total capacity possible into a network containing 
only ground stations i and./, i = 1, 2,j - i + 1, • • • , 3, and C^' (i) is the 
total capacity possible into a network containing only ground station 
i, i = 1, 2, 3. The intersection of regions (78) through (84) represent the 
capacity region upper bound for the sources of the broadcast satellite. 
We note that, if a particular ground station is included in the network 
for the purposes of calculating total capacity, then all sources com- 
municating to that ground station must be included to establish the 
region. Generalization to any G is straightforward. 

VI. SUMMARY AND CONCLUSION 

In this paper, we have reviewed various concepts proposed to date 
for introducing high-capacity, spot-beam, frequency-reuse techniques 
for communication satellites. These proposals involve varying degrees 
of complexity in an attempt to achieve a uniform grade of service over 
a service region much wider than the beamwidth of a single spot beam 
which may offer nonuniform levels of traffic. The fundamental ques- 
tion unanswered by all these approaches is the ultimate capacity limit 
of a satellite constrained in power, bandwidth, and antenna aperture 
size. 

Because the satellite system integrates multiple users, the ultimate 
capacity limit is expressed as a region over the space of transmission 
rates #i , • • • , #c into the G-ground terminals. Any particular set of 
rates contained in this region can be transmitted with arbitrarily high 
bit error-rate performance. 

Outer and inner bounds for the capacity region were derived, and it 
was shown that the shape of the capacity region is dependent upon 
both the physical locations of the ground stations and the available 
carrier-to-noise ratio. Because of the Fourier transform relationship 
between the satellite aperture illumination and the far-field radiation 
pattern, the separation between ground stations, expressed as a frac- 
tion of the narrowest beam width which can be produced by the size- 
limited satellite aperture, affects the achievable degree of frequency 
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reuse and therefore the capacity region. In the above, available carrier- 
to-noise ratio is defined to be the resulting cnr at a given ground 
station if all the radiated power available at the satellite is concentrated 
into the narrowest beam with maximum gain in the direction of that 
ground station. 

To find an outer bound, a number of independent channels or modes 
equal to the number of ground stations were formed such that the 
components of a given mode as received at the ground stations were 
orthogonal to the components of all other modes. The various modes 
exhibit differing degrees of efficiency in directing power toward the 
ground stations. The total power of the satellite was then allocated to 
the modes in accordance with a classical water-filling problem; the 
modes exhibiting better directivity are given more power. The outer 
bound so obtained is all-inclusive, that is, no tdma, fdma, spread 
spectrum, antenna pattern shaping, or scanning spot beam approach 
can achieve a greater capacity region. 

By definition, any scheme not requiring cooperating ground stations 
and for which the capacity region can be found provides an inner 
bound to the attainable capacity region since, by construction, that 
capacity region is actually achieved. The inner bounding schemes 
proposed included (i) time-sharing several maximum power-maximum 
directivity beams, (ii) formation of G independent but interfering 
beams, each with maximum possible gain in the direction of one 
ground station, and (Hi) formation of G independent beams with the 
property that each beam exhibits maximum gain in the direction of 
one ground station subject to the constraint that no interference be 
produced at any other ground station. For a given cnr, the technique 
exhibiting the tightest bound was found to be dependent upon the rate 
set ^i, • • •, 0tc that is, no one technique was found to provide the 
tightest bound over the entire capacity region. Thus, for these three 
schemes, the tightest inner bound is a composite of all three. 

For uniformly spaced ground stations, it was shown that the theory 
of discrete prolate spheriodal sequences (dpss) can be applied to 
facilitate calculation of the capacity region outer bound. The bounds 
were then applied to study, in detail, the effectiveness of frequency 
reuse for two closely spaced ground stations. A typical result shown 
there is that when the separation between the two ground stations is 
as small as 50 percent of the — 3-dB beam width of the narrowest beam 
possible, it is desirable to employ frequency reuse via noninterfering 
beams. 

Thus the capacity region bounds derived in this paper provide a 
yardstick against which the performance of any multiple spot-beam 
frequency-reuse technique can be compared. Research into antenna 
excitation techniques and multiple-user coding theory may reveal new 
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concepts yielding inner bounds closer to the ultimate capacity region 
than those discussed here. Also, the need for new devices such as 
efficient linear power amplifiers, modulators, and power combiners is 
indicated to provide higher performance, ultimately approaching the 
capacity region outer bound. 

Using the same approach as is presented here, it is possible to find 
outer and inner bounds applicable for satellites subject to more restric- 
tive antenna constraints such as phased arrays. Of course, the capacity 
region so obtained must be smaller than that presented here, which 
was derived under the more general constraint that all aperture illu- 
mination functions vanishing outside the physical dimensions of the 
aperture are permitted. 

APPENDIX 

Discrete Prolate Spheroidal Sequences 

Consider a sequence of real numbers a„, n = •••, —2, — 1, 0, 1, 
2, • • • , where ]> a'i is finite. We define the spectrum of these numbers 

as 

<*(«) = 2 a n e Jun . (85) 



Clearly, 



Moreover, 



4" + 2tt) = I a n e ilu+2 " )n = <*(«). (86) 



a n - i I *«)«*" da). (87) 

lit J_ 

We say that the sequence {a n } is bandlimited to «„ < tt if 

£(u) = 0, wo < | u | < 77-. (88) 

The energy in the sequence {a„ } is defined to be 

E = lai (89) 

n 

For a given index limit N and bandwidth w„, the dpss of order 1 is 
defined to be that sequence {a n } for which 

N 

X l (N, Uo )=-^- (90) 

is maximized, i.e., for a given index limit N, the dpss of order 1 is that 
sequence bandlimited to co which has maximum energy concentration 
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over the index range [1, N]. Moreover, for index limit N, the dpss of 
order k < N is defined to be that sequence, bandlimited to u> and 
orthogonal to each of the dpss of order 1, • • • k — 1, such that 

N 

\k(N,u„)=AJ- (91) 

2 al 

n=-oo 

is maximized. 

The energy ratios \k(N, w ), k — 1, •• •, N are eigenvalues of the 
equation 

N 

2 y(n - m)a.j, m = Aya,> (92) 

m=l 

where 

. sin u n 
y(n) = (93) 

7771 

The index-limited dpss are the eigenvectors of (92), i.e., the index- 
limited dpss of order j is 

Qj= [ay.i, •••, aj, N l (94) 

If we bandlimit the spectrum of each of (94) to to and apply the 
transform operation (87), we obtain the dpss. The dpss are doubly 
orthogonal: 

(cij > ak) = X «/." a *.n = Sjjt (95) 



X a^akjx = XjSjjt. (96) 

n=l 

Consider now the sequence b = [b it • • •, 6a/]. We wish to extend b 
outside the index range [1, N] such that the extended sequence is 
bandlimited to co and has maximum energy concentration in [1, N]. 
The result is: 23 

N 

b = £ cjOjiN, w ), (97) 

7-1 

where aj(N, w„) is they'th order dpss and 

N 

Cj = X bkdj* (98) 

*=i 
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